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Abstract. We have extended Monte Carlo simulations of hopping transport in 
completely disordered 2D conductors to the process of external charge relaxation. 
In this situation, a conductor of area L x W shunts an external capacitor C with 
initial charge Qi. At low temperatures, the charge relaxation process stops at 
some "residual" charge value corresponding to the effective threshold of the Coulomb 
blockade of hopping. We have calculated the r.m.s. value Qr of the residual charge 
for a statistical ensemble of capacitor-shunting conductors with random distribution 
of localized sites in space and energy and random Qi, as a function of macroscopic 
parameters of the system. Rather unexpectedly, Qr has turned out to depend only 
on some parameter combination: Xq = LWvoe 2 /C for negligible Coulomb interaction 
and X x = LWk 2 /C 2 for substantial interaction. (Here vq is the seed density of 
localized states, while k is the dielectric constant.) For sufficiently large conductors, 
both functions Qn/e = F(X) follow the power law F(X) = DX- fj , but with different 
exponents: f3 = 0.41 ± 0.01 for negligible and (3 = 0.28 ± 0.01 for significant Coulomb 
interaction. We have been able to derive this law analytically for the former (most 
practical) case, and also explain the scaling (but not the exact value of the exponent) 
for the latter case. In conclusion, we discuss possible applications of the sub-electron 
charge transfer for "grounding" random background charge in single-electron devices. 



PACS numbers: 72.20.Ee, 73.23.Hk, 73.40.RW 
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1. Introduction 

Electron transport via inelastic hops between localized states in disordered conductors 
has been studied for many years, with the main focus on the average transport 
characteristics (e.g., dc current dependence on temperature and applied electric field) 
and to a lesser extent on the 1/f noise - see Refs. PJ El El H] for comprehensive reviews of 
this work. The relatively recent observation 00 U\ that hopping transport may provide 
quasi-continuous ( "sub-electron" ) charge transfer gave a motivation for the extension of 
this work to the statistics of the electric charge Q carried over by the hopping current. 

The idea of the quasi- continuous charge transfer is quite simple: due to the 
electrostatic polarization, each electron hop between two localized sites inside the 
conductor leads to a step-like increase of the "external charge" Q(t), which may be 
defined as the time integral of current I(t) flowing through the wires connecting the 
conductor's electrodes to the electric field source. If an electron is transferred through 
the whole sample in one hop (as happens in the usual tunnel junctions), the charge 
step \AQ\ is equal to the fundamental charge e. However, if an electron in an extended 
conductor hops between two sites which are separated by a distance Ar much less than 
the conductor length L, then the step |AQ| is of the order of e x (|Ar|/L) <C e. (The 
exact expression depends on the sample and electrode geometry.) This means that the 
charge transport becomes nearly continuous, just as in long diffusive conductors [HE]- 
This phenomenon may have several useful applications in single-electronics, especially 
since the hopping conductors (in contrast to their diffusive counterparts) may provide 
the necessary high values of resistance R ^> h/e 2 jH] without adding too much stray 
capacitance to that of single-electron islands. 

One of the manifestations of the quasi-continuous charge transport is the 
suppression of the shot noise [31 QUI HT] . Namely, for sufficiently small values of 
the observation frequency / (with a possible exception for the 1/f noise at very low 
frequencies) the current noise spectral density Si (/) becomes approximately L C J L <C 1 
times the Schottky value 2el, where L c is some characteristic length scale. This 
prediction j^j has been confirmed in several recent experimental [T2*l ITSj and theoretical 
[HI US HH1 Ej studies of hopping. 

The goal of this work has been to study another manifestation of the quasi- 
continuous charge transfer at hopping, which is more closely related to its most 
important potential application: the ability to "ground" sub-electron amounts of electric 
charge j^j. For this, we have analyzed the simple system shown in Fig. ^ a hopping 
conductor shunts an external capacitance C with an initial charge Qj. The capacitance 
charge Q leads to a nonvanishing electric field E = V/L = Q/CL applied to the 
conductor, which causes electrons to hop through the conductor. These hops result 
in the gradual reduction of the charge Q and hence the field E. At the perfectly 
continuous ("Ohmic") conduction the process would continue until Q and E vanished 
completely (at T — > 0); however, for hopping conductors of a finite size LxW the charge 
relaxation stops at a certain finite residual charge which generally depends not only on 
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Figure 1. The system under analysis (schematically). 



the macroscopic parameters of the system, but also on the particular distribution of the 
localized sites over space and energy and on the initial charge Qi. 

Though qualitative experimental evidence of sub-electron charge relaxation has 
been obtained long ago [T%1 IT§] . to the best of our knowledge this phenomenon has 
never been studied in detail. The objective of this work has been to study the dynamics 
of this charge relaxation process, and the statistics of the residual charge theoretically 
The problem is essentially classical, but multi-particle, highly nonlinear, and statistical, 
so that most results have to be obtained by numerical (Monte Carlo) simulation using 
modern supercomputer facilities (see the Acknowledgments section below). 



2. Model 



For the hopping conductor, we have used the same model whose average transport 
characteristics and current noise had been extensively explored recently for the case of 
fixed, constant applied field E (THllIZj. Briefly, the conductor is "fully frustrated" in the 
sense that the localized sites are randomly and uniformly distributed, with a constant 
"seed" density of states z/ , over both the rectangular 2D sample of area L x W and a 
broad interval of "seed" energies e^°\ The full energy U of the system is the sum of 
the "seed" energies of all occupied sites and the electrostatic energy of the Coulomb 
interaction of the hopping electrons with each other and the external capacitance: 

v = E + i E (% - 5) - 1) « <i. *) + g- ID 

Here nj (equal to either or 1) is the occupation number of the localized site, while 
k is the dielectric constant of the insulating environment [25] • For the simplest geometry 
of a 2D conductor connecting two semi-space-shaped electrodes, the Green's function G 
in Eq. may be simply expressed as a sum over the infinite set of image charges in 
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the electrodes: 



G ( r i> r fc) = Yl 



_ \J (2nL + x k - Xjf + (y k - yjY 



(2) 



'(2nL + x k + Xjf + (y k - 

For this geometrical model, the total charge Q of the capacitor (including the 
polarization component) is 



Q = Qi- 



(3) 



where Qi is the initial charge and Xj is the j site position along the sample length L, 
while N e is the total number of electrons that have passed through the conductor from 
the start of the relaxation process until the given moment. In the limit of large charge 
(| Q | 3> Qr) the effect of capacitance on hopping transport is equivalent to that of the 
electric field E = Q/CL. 

Electron hops are permitted from any occupied site j to any unoccupied site k with 
the rate 



Ijk = r jfc exp ( — - 



where a is the localization length, and 



hT jk (AU jk ) = g- 



AU. 



(4) 



(5) 



1 - exp (-AU jk /k B T) ' 
Here AUj k is the difference of the total system energy U before and after the hop, and g 
is a small dimensionless parameter which affects only the scale of hopping conductivity 
(T = ge 2 /h. The numerical study has been carried out using the classical Monte Carlo 
technique by Bortz, Kalos and Leibowitz [21] in the form suggested by Bakhvalov et al. 
[22] , which has become the de facto standard for simulations of single-electron tunneling 
[2B] • An important feature of this algorithm is that it is not slowed down by the gradual 
reduction of hopping rates at charge relaxation. 



3. Charge Relaxation Dynamics 

Figure [2] shows, by thin lines, typical results of our Monte Carlo simulations for two 
values of the dimensionless parameter of the Coulomb interaction strength, x = e 2 z/ a/K. 
Note the logarithmic time scale and the linear scale of Q; in such coordinates the 
exponential relaxation of average charge in an RC circuit with a linear Ohmic resistor 
looks like a sharp step down at t « RC. We indeed observe such behavior at hopping 
when the initial electric field is low, i.e. in the high temperature limit. However, 
motivated by prospects of practical applications [5], our main focus is on the opposite, 
"high-field" (low-temperature) limit. Figure [2] shows that in this case the dynamics of 
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discharge through the hopping conductor is rather different: it slows down dramatically 
at Q — > 0. This is exactly what should be expected from the previous studies of variable- 
range hopping at constant applied field, which show that the hopping conductance drops 
exponentially as the field decreases [U 121 El UH UH UZj • A qualitatively similar dynamics 
is also typical for the qualitatively close (but quantitatively different) problem of intrinsic 
relaxation in electron glasses - see, e.g., recent publications [2U 1213 I2H| and prior work 
cited therein. 

It has turned out that most of the relaxation process, while the charge is sufficiently 
large (\Q\ 3> Qr), may be well described by the mean-field equation 

^ = -/ (T, E, X ) = -a (T, E, X ) EW, (6) 

where a (T, E, x) is the nonlinear conductance in the constant applied field E. In 
the low-temperature limit (fc^T -C eEr, where r is the average length of the 
hops contributing substantially into the current), we can use the following analytical 
expressions obtained by fitting the results of our numerical simulations of constant-field 
hopping within the same model [To^ I17j: 

(i) If Coulomb interaction is negligible, x 3 ^ E/E , 

1/31 



— « A (E, 0) exp 



B(£,0)§ 



(7) 



where eE a = 1/V a 2 , while A (E, x) and B(E,x) are dimensionless, weak functions 
of the applied field E and Coulomb interaction strength x- In a prior study [16], we 
have found the best fit for the pre-exponential (model-specific) function to be A (E, 0) = 
(9.2 ± 0.6) (E / E ) {0 - 80±0 - 02 \ with B treated as a constant: B (E,0) = 0.65 ± 0.02. 
(ii) If Coulomb effects are substantial, then 



— ~ A (E, x) exp 



T7 \ V2 



(8) 



For the particular value of x — 0.5, a similar approach to fitting gives A (E, 0.5) = 
(3.0 ± 0.4) (E/E ) (0 - 72±om) with B (E, 0.5) = 1.68 ± 0.07. 

For relatively low fields, E -C E Q , these formulas describe the so-called "high-field" 
variable range hopping (23 EHl E01 EI] ■ 

Broad gray curves in Fig.Elshow the results of integration of the mean field equation 
using these formulas for one value of capacitance C/Cq = 100. (The middle curves 
correspond to the best fit values, while the outer curves reflect the fitting uncertainties 
specified above.) One can see that at \Q\ ^> Qr the relaxation results may be well 
described by the mean-field approach. However, this approach does not work at Q — > 
where it predicts the complete relaxation of charge, while in reality (and numerical 
experiment) the process stalls at a certain "residual" charge. 
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Figure 2. Capacitance charge Q relaxation (at T = 0) for the cases of (a) negligible 
(x = 0) and (b) substantial (\ = 0.5) Coulomb interaction of hopping electrons. Thin 
lines show Monte Carlo results (for 6 realizations of each case) for several values of 
external capacitance C, with fixed conductor size L x W = 80 x 40a 2 . The thick gray 
curves correspond to the results of the solution of Eq. © with Eq. Q for panel (a) 
and Eq. (JSJ for panel (b) for C/Cq = 100, with the central curve corresponding to the 
best-fit parameters A and B and the outer curves corresponding to the uncertainty in 
these parameters. (See the text.) Time is measured in units of to = hv^a? / g, while 
capacitance is expressed in units of Co = e 2 voa 2 . 
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Figure 3. The r.m.s. value Qr of the residual charge at T = for negligible (x = 0) 
and finite (x = O.f and 0.5) Coulomb interaction, as a function of the conductor area 
(L x W) for different external capacitances C, and two different aspect ratios (L : W 
= 2:1 and 1:1). Each point represents data averaged over a large number (10 3 ) of 
conductor samples with vertical error bars corresponding to the uncertainty of such 
averaging. (Error bars are shown on figure, unless smaller than the symbol size). Thin 
lines are only guides for the eye. The bold horizontal line corresponds to Eq. (JHJl, while 
the bold tilted lines are the best power-law fits for large-sample data. 



4. Residual Charge Statistics 

Figure El shows some of our results for the r.m.s. value Qr of the residual charge, 
obtained for a broad range of "macroscopic" parameters of the system, including external 
capacitance C and normalized Coulomb interaction strength function of the 

conductor area L x W. (These results do not change noticeably if the systems are 
annealed after the relaxation.) 

For sufficiently small samples, the number of localized sites is so low that no 
internal hopping events may occur within the energy interval of interest, and the initial 
charge can only relax by direct tunneling between the electrodes, giving changes of 
Q in multiples of e. The Coulomb blockade theory (see, e.g., Ref. shows that at 
low temperatures such tunneling is blocked at \Q\ < e/2. If the initial charge Qi is 
random (as has been accepted in our calculations), then the residual charge is uniformly 
distributed within the range from —e/2 to +e/2, and the r.m.s. residual charge is 
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Figure 4. The same results for Qr as in Fig. El re-plotted to emphasize their universal 
scaling with system parameters. Solid lines show the best fits to the asymptotic 
behavior of Qr for large samples. 



in a good accordance with the simulation results (Fig. |3J) . 

On the other hand, if the conductor area is increased, Qr decreases, since there are 
more and more internal localized sites available for further charge relaxation. Our results 
(Fig. El) show also that Qr always increases with capacitance C and, at substantial 
Coulomb interaction, with its strength \. 

Rather unexpectedly, we have found that for a broad range of system parameters, 
all these dependencies may be very well approximated by "universal" laws, different for 
the cases when Coulomb interaction is negligible (x 3 «C Qr/CLEq) or substantial - see 
Fig. |U In the former case, Qr)c = F (X ), where 



LW C 



a 2 C 

while in the latter case Qr/c 
LW 



= F X {X X ), 
LWk 2 



where 



x 2 c 2 



c 2 



(10) 



(11) 



At small values of their arguments, both functions F tend to l/vT2, in agreement with 
Eq. (J?|J). Their asymptotic behavior is also functionally similar, F(X) — > DX~^ at 
X — > oo, but with different best-fit values of the coefficients: for \ = 0, D = 0.64 ±0.01 
and (3 = 0.41 ± 0.01, while for x ~ 1, D = 1.1 ± 0.1 and (3 = 0.28, with the error 
bar about 0.03 for the dependence on C and of the order of 0.01 for other variables 
contributing to X x . 
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5. Discussion 

For the case of negligible Coulomb interaction, the asymptotic power law for function 
F (X ) may be readily explained , using the basic ideas of the Coulomb blockade 0. 
Charge relaxation continues with the reduction of the system energy (on the average, 
dominated by the capacitor energy U) until the number N of localized sites available for 
hopping becomes less than one. If the capacitance charge before a hop is Q, the range of 
capacitive energy of available initial sites is AU ~ Q 2 /2C, so that the average number 
of such sites per unit area is ~ u AU ~ u Q 2 /2C, and their total number in the 
sample of area L x W is Ni ~ LWni ~ LWvqQ 2 /2C ' . In order to estimate N, we need 
to multiply iVj by the average number Nf of available final sites for each initial site. For 
small changes of charge, \AQ\ <C e, the area |Ax| x W where such states can reside is 
much smaller than the sample area L x W, because such charge change corresponds to 
a hop by distance \ Ax\ = Lx\AQ\/e < L. Hence N f ~ LWis (\AQ \/e)(Q - AQ) 2 /2C 
and we get the following estimate 

N ^ N , Nf ^(^) 2 9^9m^9L. (12) 

V I e 

Now, from the natural requirement that N drops below 1 as soon as \Q\, \AQ\, and 
| Q — AQ\ all become, on the average, of the order of Qr, we get 

^(^£!)- 2/ W's, (13) 

which when compared to the power law F(X) discussed above gives (3 = 2/5 = 0.40, 
i.e. inside the narrow interval 0.41 ± 0.01 given by the numerical experiment. 

For the case of substantial Coulomb interaction of hopping electrons, the situation is 
more complex - see, e.g., the discussion on pp. 435-443 of Ref. PJ. It is well documented 
that "external" transport (bringing electrons into and out of the hopping conductor) 
may be well understood in terms of the simple quasiparticles introduced by Efros and 
Shklovskii [2], with energy 



e 3 =e 



K 



In 2D systems, their density of states at low energies e is given by the famous Coulomb- 
gap expression [2j 

u(e)^—\e\. (15) 

If we naively repeat the above calculation of Qr, just replacing u with u(e) from the 
last expression, we get 

Qr [LWk^ - 2/g 
e \ C 2 

i.e. the experimentally observed universality (X x = LWk 2 /C 2 ), but with an exponent 
(3 = 2/9 ~ 0.22 which is significantly outside of the experimental interval 0.28 ± 0.01. 



X x ~ 2 '\ (16) 
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Actually, for intra-sample transport, more adequate quasiparticles may be the 
so-called "dipole excitations" (essentially, electron-hole pairs with correlated energies) 
whose density F(Q, r) depends on both the pair energy Q and the distance r between 
the pair components (see [3] p. 435). In contrast to constant-field transport, the residual 
charge statistics are dominated by large-size pairs (hops), with a;— component of the 
order of L(\AQ\/e) and y— component of the order of W. If we neglect, for such hops, 
the interaction of the pair components in comparison with Q, then F depends only on 
energy: 



For energies f2 much less than both the cutoff energy A and the Coulomb gap width, 
this integral yields 



Now, following the arguments used above, we can accept Q ~ Q 2 /2C, and take LW 
for the possible area of the pair centers, and L(\AQ\/e)W for the pair area. After the 
integration of F from to Q, for the possible number of pairs within our energy range 
we get 



Again, requiring that iV ~ 1 at Q, \AQ\ ~ Q R , we get back to the estimate given by 
Eq. ffToT). 



It is not quite clear presently whether the discrepancy between these analytical 
arguments and the results of our numerical experiments may be overcome by an account 
of electron-hole pairs of smaller size, with strongly interacting pair components. 

6. Offset Charge Grounding 

The results of this work allow one to estimate the prospects of applying hopping 
conductors as "grounding" devices for the random background charge in single-electron 
devices. Figure El shows this idea on the example of a single-electron transistor 
[ElE]. Charged impurities, randomly located in the vicinity of the transistor's single- 
electron island, induce on it a net polarization charge. The "integer" (e-multiple) 
part of this "background" charge is automatically compensated by tunneling through 
the transistor's tunnel junctions, but its fractional part — e/2 < Q < +e/2 cannot 
be compensated in this way. This random charge is equivalent to a random shift 
AV g = Qo/Cg of the gate voltage; such shifts are one of the main obstacles on the way 
toward integrated circuits using single-electron devices, because for most of them the 
tolerable background charge range is as narrow as ~ O.le [Oj. The problem may be solved 
by connecting the single-electron island to "ground" through a hopping conductor which 
would provide a slow relaxation of the background charge jH] . (For digital applications, 




(17) 




(18) 




(19) 
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Figure 5. Background charge "grounding" using a hopping conductor (schematically). 



the characteristic relaxation time has to be much longer than at least the circuit clock 
cycle, and more preferably the full time of the calculation performed by the circuit.) 

For typical hopping conductors technologically compatible with silicon technology 
(e.g., amorphous semiconductors and metal oxides), the dielectric constant k is of the 
order of 10, while the electron effective mass m ~ 0.2m . This gives the localization 
radius a ~ h 2 n/me 2 ~ 3 nm and the level splitting scale e 2 /na pa me A /K 2 h 2 ~ 30 meV. 
In order to stay on the dielectric side of the metal-insulator transition, the average 
distance between the localized sites should be above ~ 4a jT]; for the 3D density of 
states u 3 this gives the condition u 3 < 10 19 eV _1 cm -3 . This condition is well satisfied, 
e.g., for most species of amorphous silicon, where z/3 at mid-bandgap is of the order 
of 10 16 eV _1 cm~ 3 (see, e.g., Ref. [32J). For thin films of such material with thickness 
t ~ a ~ 3 nm, the 2D density of states z/ ~ 3 x 10 9 eV _1 cm -2 . For these parameters, 
the Coulomb interaction parameter \ is much smaller than 1, and we can use Eq. (jTBj) 
for estimates. Even for the least demanding applications of single-electron devices, 
the electron addition energy e 2 /C should be at least 30 &bT jH], so that according to 
Eq. (fT3*j). X has to be above ~ 300. 

Let us accept L = W in order to minimize the conductors' self- ("stray") 
capacitance C s (which, as we will show shortly, may present a major problem) at fixed 
area L x W. For the usual conditions of low-temperature experiments with single- 
electron devices, T ~ 0.1 K, C may be of the order of 10~ 14 F, so that with our 
parameters L should be above ~ 30 nm. This is less than the typical length (~ 1 
/im) of the single-electron island in such experiments, so that the grounding idea may 
actually work [3*3*] . 

On the other hand, for the most important case of room-temperature single-electron 
devices (T pa 300 K), the island capacitance should be much less, C < 10~ 18 F, so that 
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the quasi-continuous conduction is only possible at L > 15 jum. Stray capacitance C s 
of such a conductor would be larger than ~ 10~ 15 F, i.e. much larger than C, thus 
increasing the total effective capacitance of the island well above the acceptable value. 

To summarize, our calculations indicate that the fractional charge grounding is 
possible, but practicable only for low-temperature experiments rather than for room- 
temperature single-electron devices. Fortunately, by now an alternative way to solve (or 
rather circumvent) the random background charge problem in digital nanoelectronics 
has been suggested. This approach is based on reconfigurable hybrid CMOS-nanodevice 
digital circuits which may be re-routed around "bad" devices - see, e.g., Ref. [33]. Recent 
calculations have shown that this approach may provide defect tolerance up to ~ 10% 
in memory circuits and > 20% in logic circuits. This is much higher than the estimated 
lower bound on the fraction (~0.1% [9J of single-electron devices whose threshold is 
substantially shifted by random background charges. 
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